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Abstract 



In this paper we develop a formalism to evaluate wave functions in momentum and 
coordinate space for the resonant states dynamically generated in a unitary coupled 
channel approach. The on shell approach for the scattering matrix, commonly used, 

(f) ■ is also obtained in Quantum Mechanics with a separable potential, which allows one 

to write wave functions in a trivial way. We develop useful relationships among the 
couplings of the dynamically generated resonances to the different channels and the 
wave functions at the origin. The formalism provides an intuitive picture of the 
resonances in the coupled channel approach, as bound states of one bound channel, 
which decays into open ones. It also provides an insight and practical rules for 

^ evaluating couplings of the resonances to external sources and how to deal with final 

state interaction in production processes. As an application of the formalism we 
evaluate the wave functions of the two A(1405) states in the ttT,, KN and other 
coupled channels. It also offers a practical way to study three body systems when 
two of them cluster into a resonance. 

1 Introduction 

The chiral unitary approach to hadron dynamics has brought a new perspective to deal 
with the interaction of hadrons and the nature of some resonant mesonic [1-7] and baryonic 
states [8-13, 15-17] which appear dynamically generated from the interactions and, thus, 
have a nature quite different to standard qq states. With some different formulations at the 
beginning [6,8,9], the more recent work uses the on shell formulation firstly established on 
the basis of the N/D method in [18], where the potential and the t-matrix in momentum 
space factorize outside the loop function implicit in the Bethe Salpeter equation in coupled 



channels that one uses in those approaches. This is a very practical way to deal with the 
problem since one renders the coupled integral equations into a set of algebraic equations, 
paying a small price which is the fine tuning of some subtraction constant appearing in 
the dispersion relations involved. The approach is practical and useful, but carries also 
a handicap which is that one deals with amplitudes in momentum space, and couplings 
of the dynamically generated resonances to the different channels, and nowhere do wave 
functions in coordinate space appear in the approach. For instance, all properties of a 
resonance are given in terms of the mass and width and its couplings to the different 
channels, obtained from the residues of the amplitudes at the poles in the complex plane. 
Any intuition about the wave function of the different channels, its magnitude and extent 
in space, is lost in the approach as well as the meaning of the resonances and the discrete 
values obtained for the resonant energies. Obviously, these are magnitudes that help to 
understand the microscopical composition of the dynamically generated states and, hence, 
a most welcome information. However, this is not all, since the scattering amplitudes do 
not contain all the information of the wave function. They reflect the wave function at 
long distances. For some observables the wave function at small distances is needed. This 
is the case when one studies the response of states to external sources, which require to 
evaluate form factors, or expectation values of different observables. The wave functions, 
whether in momentum space or coordinate space, are then needed. To study couplings of 
the resonance to local sources the wave function around the origin is required, but in the 
study of form factors one needs to know the wave function at all distances. We provide 
them both in the present work. They are already proving very useful to study new systems 
with three or more particles, when two of them cluster to form one of these dynamically 
generated resonances [20]. 

The first steps in this direction were done in [19] in order to understand the X(3872) 
resonance in terms of a molecule of D°D* and D + D*~ and their charge conjugates. The 
reason is that the D°D* component is slightly bound, while the D + D*~ one is bound by 
about 7 MeV. This has as a consequence that the D°D* component stretches up to very 
large distances, while the D + D*~ one is more confined, and this has important repercussion 
in the interpretation of observables regarding this resonance. 

The work of [19] provided also an extension of the two channel (D°D* and D + D*~) 
formalism to a general one with many coupled channels, but all of them bound. Hence, 
in [19] one deals with bound states which as a consequence have a discrete spectrum of 
energies, and the wave function was only evaluated for the eigenstates of the system. It 
is most advisable to extend the approach to the case where one has also coupled channels 
which are open at a certain energy, which is the majority of the cases in the studies done 
with the chiral unitary approach. One might think that the problem is formally identical 
to the one of the bound coupled channels, but, although there are certainly analogies, there 
are also subtle, and important, differences which call for a detailed study. The first big 
difference is that now, one has wave functions, those of the open channels, that extend to 
infinity, while the wave functions of the bound channels will remain constrained in space. 
The wave functions for the open channels will not be finite normalizable, while those of the 



bound channels with remain finite normalizable. This problem has been found before and 
the difficulties of dealing together with bound states and open, unnormalizable states, has 
been pointed out [21,22]. We have found a way to deal with this problem and answered 
questions like: What are the magnitudes that matter in this case when it comes to evaluate 
observables in terms of the wave functions? What is the meaning of the resonances in this 
approach, and why the number of resonant energies are finite? These are some of the novel 
questions that one may now ask and which have to be properly addressed. This is the 
purpose of the present work. We shall develop the formalism and then we will apply the 
results to evaluate the wave functions in the ten coupled channels of the two A(1405) states 
obtained in [15,23]. 

2 Formalism 

We follow closely the formalism of [19] adapting it to the case of open channels. First 
we will consider the relation between a coupling constant and the wave function in one 
channel case. In a second step we will extend the formalism to the case of multiple coupled 
channels. 

2.1 One channel case 

As in [19], we take a potential V in s-wave as a separable function in momentum space 
with the modulating factor being a step function, 

(p'\V\p) = V(p',p) = v6(A-p)6(A-p) (1) 

where p and p' indicate the moduli of the three momenta p and p', A is a cutoff in 
momentum space and 9(A — p) is the step function. In Eq. ([[]), v does not depend on 
momentum. It is worth stressing that the use of this simple form is sufficient to develop 
the formalism and it has the advantage that it leads to the same on shell factorized Bethe 
Salpeter equations (here we shall use their non relativistic Lippmann Schwinger form) than 
in the chiral unitary approach. 

2.2 The wave function 

The Schrodinger equation reads 

{H + V)\i>) = E\iP) (2) 

where H is the kinetic term, V is the potential, ip the exact wave function of the full 
Hamiltonian H = H + V and E the energy. Since we have open channels we have solutions 
for any value of E. In order to derive the Lippmann Schwinger equation we proceed as 
usual, introducing the solution <p of the kinetic energy Hamiltonian H for the same energy 
E. 

H \(f>) = E\<f>) . (3) 



From these equations, we obtain 



(E-H )\^-4>) = V\il>) (4) 
1^-0) = ^^H^W) ( 5 ) 
' W) (6) 



E-H, 



and the wave function in momentum space, 

(pm = w>+ / d 3 p'd 3 p"{p\ ¥ ^ w y){p'\v\p")( P "\ 



d{^-p) .. / ,3 /////, 



= W)+ F m 1, rfi/2„4.if V WWW) ( ? ) 

£/ — mi — m 2 — p^/Zfi + xe J p//< A 

where \i is the reduced mass of the two particles. In the bound state case (E < mi + m 2 ) 
the term (p\4>) does not appear, E — mi — m 2 — p 2 /2fi cannot be zero for any energy and 
for some energies one finds discrete eigenstates. However now, any energy E > vrt\ + m 2 is 
allowed. Let us integrate Eq. (JTj) over p, 



/ d 3 p(p\ip) = 6 (A - sj2^i(E -m x - m 2 ) 

Jp<A ' ^ 



J P <A E-m 1 -m 2 - P 2 l1[i Jp» <A 



(A - y/2(i(E ~m 1 - m 2 )) + Gv [ d 3 p"(p 



where 



G= I d 3 p- - . (9) 

J p<A E - mi - m 2 - p z /2fi + le 

the term 9 ( A — ■s/2/j,(E — mi — m 2 ) ) in Eq. (J8]) comes because of the normalization that 

we impose on the states 

(p\p') = 6 3 (p-p') (10) 

since \<f>) = \p') such that p' /2/i + mi + m 2 = E. Since the integration in Eq. ([9]) 
has the singularity at E = mi + m 2 + p 2 /2/i, we put +ie which guarantees an outgoing 
solution for the Lippmann Schwinger equation. In what follows, we will assume A > 
\ l /2[i(E — mi — m 2 ). This is essential to get the unitarity properties from the G function 
of Eq. ([9]) and we must be certain that this occurs for all channels in the coupled channel 
case. From Eq. (jSJ) we have now 

(1 - Gv) f d 3 p(p\iP) = 1 (11) 

Jp<A 

d 3 p(p\i>) = T^7. (12) 



P<A 



1-Gv 



and if G is complex this means that f A d 3 p(p\ip) is complex and hence ip\ip) is necessarily 
complex. 

From Eq. ([7]), the wave function in coordinate space can be equally evaluated by 

(a#) = d 3 p(x\p}(p\ilj} 

— e ip'x 



(2vr) 3 /2 

[ 1 1 1 

+ J p<A d3p j2^ 6tP * E-m 1 -m 2 -pi/2 l i + i^Y^Gv (13) 

and we can obtain the value of the wave function at the origin in coordinate space 

Now we define if) = (27r) 3//2 ?/>(0) and we obtain 

$ = (2tt) 3 />(0) = r= ^- . (15) 



and by means of Eq. ( 1T2|) 

$ = / d 3 p(f#) (16) 

Eq. (JTBjl is also found in the case of bound states, but Eq. ( fl5|) is new for scattering. For the 
case of a bound state one finds in [19] ip = gG, where g is the coupling of the bound state to 
the channel considered. This difference must be clarified. In Eq. ([3]), we considered a plane 
wave (J), which is needed to satisfy the boundary condition at infinity. For the bound wave 
function, <f> does not appear because (E — H )\<p) = has no solution for E < mi + m<i. In 
this case the Schrodinger equation has only a few discrete eigenenergies E a . 

2.2.1 The coupling 

Let us consider the Lippmann Schwinger equation. Let us define T such that T\4>) = V\ip). 
From Eq. (EJ), we obtain 

t = v + v ^zrjr o T ■ ( 17 ) 

From Eq. (fT7|) . we can write 

(p\T\p'} = (p\V\p>) + [dY- ^^^ „ 2/0 (P"\T\P') (18) 

J E — rrii — m 2 — P /2/i 
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which has the solution 
and then 



(p\T\p') = 0(A - p')6(A - p)t 

v 1 



t = v + vGt 



(19) 



(20) 



1-vG v- l -G 

For real E we do not have poles since vG is complex for real E. However, when we have 
1 — Re vG ~ we have a near pole and hence an enhancement of t, and thus a resonance. 
As we know, we can even find a pole in the complex plane in the second Riemann sheet 
associated to this resonance. 

In the vicinity of the pole, t also can be written as 



t 



9' 



1 



E - E R + iT/2 



G 



(21) 



We define the resonance energy E R when one satisfies v l — Re G(E R ) = 0. Using the 
equation below, and assuming v constant as a function of energy 



v- 1 - Re G{E) 
we rewrite Eq. (|2"T|) as 



-1 



Re G(E R ) 



dReG 



BE 



{E - E R ) + 



(22) 



(9Re G 



BE I-&R 



ER {E-E R )-iImG{Ei 



(23) 



which gives us the two conditions 
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r 

2 



~\ dE )e r 
-g 2 lmG . 



(24) 
(25) 



The integral for the G function defined in Eq. (j9]) can be performed analytically and 
we obtain 



G(E) 
k 



k A — k 
-87r/i(A + - In -) - i47r 2 nk 



(26) 

(27) 



\r2JJ(E -mi- m 2 ) ■ 

It is instructive to check that an analytical extrapolation below threshold of Eq. fT26l) . 
putting k = ±«7 (7 = y/2fiB, with B the binding energy), leads to the formula for the G 
function for bound states (Eq. (27) of [19]). The above equation gives us the coupling g as 

(28) 



i/J 2 TC 



In (A - fc R )/(A + k n ) 1/2 - A; R A/(A2 - fcg 



The first thing we observe is that g 2 < and T < 0, which is not physical. The result 
obtained simply expresses the fact that a constant potential v as a function of energy does 
not lead to a resonance (an energy dependent v could give a resonance). However, with v 
negative one can obtain bound states as found in [19]. A resonance state in one channel is 
usually associated to a barrier in coordinate space which is not reproduced by v constant. 
The situation will be very different when we go to coupled channels. 

The other feature we must to mention is that the limit of g 2 when k-& — > is given by 

^ = "«£% (29) 

which is not zero, unlike the case of bound states [19], or the case that we shall see in 
coupled channels, where the existence of only one bound channel guarantees c? 2 — ^ for all 
channels including those with energies above the threshold of the bound channel. 

We can also obtain the asymptotic behavior of the wave function at r — y oo from 
Eq. (1T31 and we find 



D ip r 

(27r) a/2 (a#> > e tp ' x - An 2 fi — "—'— (30) 

1 vG r 



which tells us that the scattering matrix f{6) of Quantum Mechanics is given by 



f{9) = -Ait l n T ^ = -ATt i nt . (31) 

2.3 Coupled Channels 

We extend now the formalism to N coupled channels, where at least one of them is bound. 
We take again 

(p'\V\p) = 9(A - p)9(A - p')v (32) 

where v is now a. N x N matrix. 

2.3.1 The wave function 

The first thing when dealing with coupled channels is that one must find the boundary 
conditions for the physical process that one is studying. If we wish to create a resonance 
from the interaction of many channels at a certain energy we must take a channel which is 
open at this energy and make the two particles collide, starting from an infinite separation 
at t = — oo. Let us call channel 1 to this open channel that undergoes the scattering. 
The equations to solve, with the appropriate boundary condition of a scattering state for 
channel 1 are 

m = 10) + ^rjf/^) ( 33 ) 



where 



' IV*> ' 



\<k) 






(34) 



and |0i) = \p'), such that p' /2/ii + Mi = E, where we will use the notation M{ = mii + m2i 
and fii = m li m 2i /(m li + m 2i ). 

Following Eq. (J7J), the wave functions in momentum space are written 



(p|Vi) = (p\(t>i) + 

(p\A) 



9{A-p) 



E-Mi- p 2 /2n 
9{K-p) 



-— E^7 <*vwi> 

1 + «e i -" Jp"<:A 



£ 



dVVU-) (<^i). 



E-Mi- P 2 / 2 /^ + ie— j r<A 
And integrating over p, we have 

/ d 3 p(Pl^i) = l + GuX)^/ d VVU'> 

Jp<A ■ Jp"<A 

[ d 3 p(p\^} = GuJ^Vij f dV(p"\^} («^1) , 

./p<A ,• Jp"<A 



(35) 
(36) 

(37) 
(38) 



where once again we assume A to be bigger that the on shell momenta of the particles for 
all the open channels. In Eqs. ( 157)) . ( 155)) the diagonal G matrix is given by 



(G 1 



G 



Go 



\ 



(39) 



V G N J 

with Gj given by Eq. (Q for each channel. Thus, we have 

V>i = 1 + G11 7, ^ij^j 

j 
where ipi = J <A d 3 p{p\tpi) and Eq. (|4ip is then written in matrix form as 



(40) 
(41) 



'f 



M = (1 - Gv)-' I 







(42) 




v / 



which provides the values of ipi as 

A = (1 " Gv)-, 1 . (43) 

The scattering matrix in coupled channels is given formally by Eq. ( TT71) . but T is now a 
N x N matrix. It has also the form of Eq. ( Tl9l) and now the N x N t matrix is given by 



t=[l- vG}- l v = [v- 1 - G}- 1 . (44) 

By means of this equation, the combination vip appearing in Eqs. (135]) . fl36|) can be written 



as 



Viji>j = Vij(l - Gv)^ = [v(l - Gv) ^i 

= [w-V 1 -GO" 1 k = [(«" 1 -G ! )" 1 ]ii = *ii (45) 



Taking this into account, the wave function equations of Eqs. ( |35l) . ( 1361) become 

<*'•> = w^ + g-^/k + fe '"'^ (46) 

w *> = E-M<-7L+« ta(E) (M1) (47) 

and we can write the wave function in coordinate space as 

Wi) = jA^^+l d3 PjAm eip % M 1 2/ 9 .. *ii(g) (48) 



( *"> = X < / P p^ e "' %-M.-p72.. + ^ " (fl) (M1) (49) 



with |p'| = ^2^ X {E-M X ). 

Once again we can make the limit r — > oo and we find the asymptotic solutions 

ip'r 

(27r) 3 / 2 (^|^i> ► e^* - 4ttV— *ii (50) 

(2ir) 3/2 (xlipi) > -4:7i 2 fii tn (open channel) (51) 

r— >oo r 

— 47r 2 /ij tji (bound state) 

r 



where fcj = y/2jXi{E — Mj) for open channels and «$ = a/2/Uj | -E — Mj | for the bound 
channels. 

As we can see, the channel 1 contains the plane wave from the scattering state while the 
other channels only have the wave functions generated from the collision of the particles 
of the scattering channel. For bound channels one generates a bound wave function while 



in the other channels we have an outgoing wave. Eqs. ( l5Tj) show the connection of the 
scattering matrix t with that of Quantum Mechanics fn, such that 

(27r) 3 / 2 (z# 4 > ► [E fil{ e) e — (52) 

which implies 

sL=£i'«<«>i"- (53) 

Thus we have 

f a {6) = -i^y/jHf^tn (54) 

with no angle dependence since we are dealing with s-waves from the beginning. 

We can go back to Eqs. ( HH1) . (j4"9l) . by taking (27r) 3 / 2 (a; = 0\ipi), which is ^ (see 
Eqs. (fl5l) . (TT6l)). and calling ?/>} ou the outgoing wave of (a?|"0i) at a; = (term of J d 3 p in 
Eq. (liHjl ). we obtain 

^ out) (£) = G^fytniE) (55) 

^(£) = Gi(E)tn(E) (56) 

which generalizes similar equations obtained for the bound case in [19], but substituting 
tn, tn by gri, gi respectively. Here E is a continuous variable and we see that the wave 
functions depend on the energy, but for resonant energies the wave functions at the origin 
grow roughly like t, since G(E) is smoother than t at the resonant peak. Thus, we get an 
intuitive idea about the meaning of a resonance, which is that for resonant energies there 
is an accumulation of strength of the wave functions at the origin. We can also see that 

_G A %_ = U,^g 1 

the last part of the equation holding if we can represent the tn amplitude approximately 
as0tfi/(V5-Af fl + ir/2). 

We based all our approach on a potential with a sharp cut off of Eqs. ([T]) and fl32l) . The 
work can be easily generalized to the use of a potential of the type 

(p'\V\p)=vf(p)f(p') (58) 

as also done in [19]. The modifications are minimal and can be followed from Ref. [19], 
section VII. The most important for our discussion is that the function ip that represented 
the wave function at the origin is replaced by 

$*> = Jd 3 kf(k)(m) 

d 3 xi>i{x)f{x) (59) 



where f(x) is the Fourier Transform of f(k). Thus, the wave function at the origin in the 
present approach is replaced by the folding of the wave function with the Fourier Transform 
of the factor f(p), in practice an average of the wave function close to the origin. 
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2.4 Width of a resonance and partial decay widths 

Within the chiral unitary approach, a different convention for the T matrix is commonly 
used. The relationship between these matrices was given in [19]. Since many relationships 
are familiar to practitioners of the chiral theoretical approach in the field theoretical nota- 
tion, we adopt in this section the latter notation and one has for the case of meson baryon 

states 

7^ T = 32TT 3 ^JI^ v ^(2M Bl 2M Bl )~ 1/2 t il 

= -8ir v G(2M m 2M Bi )- 1 / 2 f tl (6) (60) 

where M B i are the masses of the baryons. 

The optical theorem in the scattering of channel 1 is stated as (k\ = p') 

" fcl ^f>i/a(0)| 2 ^ (61) 



M Bl f-f "" v n k 

2=1 



where we have made use of Eq. ( |53p . the Air comes from the dfl integration and N' are the 
numbers of open channels. By using Eqs. ( |54|) . (1601) we find 



A" 



Im ^=-^ET |T " |!f5 • < 62) 

1=1 * 

Now assume we have a resonance in the coupled channel problem and hence close to 
y/s = Mr we have 



2 



FT _ 9\ . -r t,FT 



7 2 



9\ 



T FT - *± ■ Tm T br - — iL- (63) 

11 ~ ^-m r + 1 t/2 ' imil1 vs=m r - r/2 ' {b6} 

On the other hand, 



t-iFT _ 9l9i . \ T FT\2 



^ 2 » a (64) 



yfi-M K + iT/2 ' ' iJ ' V~s=m r (r/2) 2 

with ^ the coupling of the resonance to channel i in the field theoretical convention. 
By means of these equations we find 

N' N' 

- = Y-^g l % = Y^ (65) 

2 ^ 4vr Mr ^ 2 K ' 

i=\ tt i=l 

the last equation holding since in this notation the partial decay width of the resonance to 
an open channel is given by 

9i h . (66) 



2vrM R 
11 



2.5 Coupling of a resonance to external channels and final state 
interaction considerations 

Sometimes one wishes to evaluate the coupling of the resonance to an external channel, 
where the coupling is sufficiently weak not to deserve to be taken into account as one of 
the coupled channels in the approach. Since the resonance has been created by coupled 
channels, the coupling of this external channel to the resonance will come from its couplings 
to the different channels. We show the simple case where the external channel couples to 
the building channels through a zero range interaction and let tj ex be the transition matrix 
from any channel to the external one. The coupling of channel 1 to the external channel 
can be depicted by the series of terms of Fig. [TJ We can ignore the first term of the series 



Figure 1: 



if we are only concerned about the resonance contribution. Hence we get 

N N 

tl „ = E tl . G ^E gf GA . (67) 

i=l i=l v ' ' 





which we would like to equate to 

tlex = _ 9l9cx (68) 

where g ex is the coupling of the resonance to this external channel. We readily obtain 

N 

9c x = ^ 9i G i f ic X ■ (69) 

i=l 

Note that we have assumed tj ex to be of zero range, hence constant in momentum space. 
Since in the series of Fig. [1] we have a v8(q — A)8(q' — A) in each four leg vertex, except 
for the last one, one guarantees that in the loop function one is implementing the cut off 
A, including the last loop. This is interesting to note because in Field Theory, in principle 
the last loop can be regularized in a different way. However, in the Quantum Mechanical 
approach we see that all loops are regularized with the same cut off and the last loop is the 
same one appearing in the scattering problem, provided the range of £; ex is shorter than 
1/A. The result of Eq. (169]) is often used in problems using the chiral unitary approach [25] 
and here we find a justification for it. 

The other subject worth discussing from the present perspective is the final state inter- 
action. Assume we have a physical process in which one state is produced. We assume that 
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the production, for instance in a weak process, is of zero range. Let us then assume that 
this state couples to iV different channels through strong interaction. On top of the direct 
production of the original state, say channel 1, one could have as well the direct produc- 
tion of any of the coupled channels, which then will make a strong interaction transition 
to channel 1. This is depicted in Fig. [2j 

<+o< i +o6< i +- 

Figure 2: 
The full production amplitude will be given by 

N 
Pl = Pl + J2 P i G ^ ■ ( 70 ) 



i=l 



Since Pi is considered of zero range, and hence a constant in momentum space, Gi appearing 
in Eq. ( ITUl) is regularized with the same cut off A in Eq. ( ITUl) than in the scattering problem. 
This is because all the four leg vertices appearing in Fig. [2] contain v9(A — q)9(A — q'), 
which impose the A cut on the loop function. In principle, in Field Theoretical studies the 
first loop to the left in Fig [2] could be regularized in a different way than in the scattering 
problem. Yet, the Quantum Mechanical treatment shown here reveals that one can use the 
same loop function as in the scattering problem provided the primary production vertex 
corresponds to a short range process, shorter than that implied by the cut off A used in 
the scattering. Eq. (ITU)) is often used in calculations using the chiral unitary approach [26]. 
In particular, if we have only one channel, Eq. ( ITU]) gives 

Pi = Pi(l + G 1 t 11 ) = Pi— , (71) 

which is consistent with Watson's theorem [27] and it was also discussed in Ref. [6]. This 
is the most popular way to implement final state interaction, simply multiply with tu 
the primary production amplitude, normalizing data at one energy, thus, using only the 
energy dependence of tn in the analysis [28]. Once again, the present approach shows the 
conditions upon which such procedure can be trusted; i) Dominance of just one channel, 
ii) short range nature of the primary production, Hi) very weak energy dependence of v 
with respect to t. This lather condition is flagrantly violated in the case of 7nr production 
around the <r(600) region, where v depends more strongly on the energy than t, in which 
case, keeping vn in the tu/vu factor is essential. This was shown in [29, 30] and it 
was essential to interpret the apparently narrow tttt "cr" structure in the J/ip — > uj , k + 'k~ 
production [31] in terms of the generally admitted wide a resonance [1,32]. 
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2.5.1 The couplings for resonant states 

The scattering matrix is given by the obvious generalization of Eq. (|20l) for coupled channels 
and it can be rewritten as 



t 



Av 



det(l - vG) 



(72) 



with A = [det(l — vG)](l — vG) 1 . Now there are open channels and any energy E is 
allowed. E is a continuous variable. Technically we can write 



Av 



Av 



det(l - vG) Re[det(l - vG)} + zlm[det(l - vG)] 



(73) 



and look for E R such that Re[det(l — vG(E R ))] = 0. Around the resonance energy E R , we 
can write 



Re[det(l - vG(E))] = Re[det(l - vG(E R ))} + — Re[det(l - vG{E))] 

oE 



E=Ei 



(E - E R ) + 



and then 



Av 



dE J 



Re[det(l - vG(E))] {E - E R ) + ilm[det(l - vG{E R ))} 



(74) 
(75) 



E=E n 



In the vicinity of the pole, t can be written as 






9i9j 



E- E R + iT/2 

and from Eqs. f!75|) and f!76|) we obtain the couplings and width 

{Av)ij 



(76) 



9i9j 

9i_ 
9.) 

r 

2 



^Re[det(l -vG(E))} 

{Av) tJ 



E=Eu 



JJ 



E — E-a 



(Av, 
Im[det(l - vG(E))] 



_d_ 
8E 



Re[det(l - vG(E))] 



e=e^ 



(77) 

(78) 
(79) 



Let us note that some equations relating the couplings obtained for bound states in [19] 
were a consequence of the dominance of just one eigenstate for the discrete eigenenergies 
in a sum over intermediate states. This does not hold here, where E is a continuous 
variable. Also other properties where based upon imposition that the total wave function 
was normalized to unity, which is no longer possible here. Yet, there is one interesting 
property which still remains. Assume we have a resonance state close to the threshold 
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of a bound channel, a. In this case dG a /dE — > — oo, then —vdG a /dE dominates in the 
denominator of Eq. ( 177|) and all the couplings go to zero, as was also the case in the bound 
states [19], and in agreement with the claims made in [33]. 

One can also see that under the dominance of one channel where \v aa \ 3> \vij\ (i, j ^ a) 
and v aa negative, one can find a resonance energy below the threshold M a , and if this is 
close to threshold, such that (dG a /dE)\ E=El{ dominates over all the other terms, one finds 
approximately, 

10) 



1 



9a 



rs_/ 



dE G a(E) 

E=E R 

like in the case of one channel bound state [19]. We will come back to this case, which offers 
a nice interpretation of the resonance. It corresponds approximately to a bound state of 
a chosen channel (hence justifying the discrete values of the resonant energies) which can 
decay into the open channels, an intuitive picture of a resonance, which corresponds very 
approximately to most of the dynamically generated states in the chiral unitary approach. 
It is also interesting to obtain the relationship of the couplings and the wave functions 
at the origin. For the case of bound states, it was particularly simple and we found_| 
[19] giGft = ipi, where a refers to G calculated for the energy of the bound state. This 
relationship is tied to the discrete spectrum of the bound state and does not hold here. 
Instead we found Eq. (H3"|) . which we can rewrite as 

^ = (1 - Gv)£ = [v-^v- 1 - G)- X ] a 

= (y-H) a = (v -1 )«-f — — t • ( 81 ) 

The last equation relates the couplings to the wave functions at the origin. The value ipi 
of Eq. ( )8T]) is tied to the choice of normalization made by us in Eq. (TTOj) . A more useful 
relationship, independent of the precise normalization, is given by 

'' ; - {V ~ l)ij9j -> ^ = a(v-% 9j (82) 



which can be equivalently cast as 



g% Vijipj 



3) 



9k V kj 1pj 

In the particular case where one channel a is dominant over the others, \v aa \ ^> |i>y 
(i, j 7^ a) and \v ai \ 3> |i)y| (i, j ^ a), Eq. (183|) reads as 



'i 



9i Vi a 1p a Vi a /oa\ 

— = r = — • ^ 4 J 

9a V aa 1p a V aa 



1 The reader will note that in Eq. (118) of [19] there are two equations that should be separated by a 
space but they are not (the some happens in Eq. (120)). 
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2.6 Intuitive picture of a resonance as a bound state of a channel 
decaying into open ones 

Let us assume that we have just one bound channel, a, and we look for bound states of 
it. We will solve the Schrodinger equation and find different discrete eigenenergies E a 
for different eigenstates ip a . Let us just take one of these states and let us assume that 
\v aa \ ^> \vij\ (i, j 7^ a). The picture is that we have essentially a bound state that decays 
into the channels i ^ a (assume the rest of channels open). An extreme picture of that 
would be an electron in an excited state of an atom which decays emitting a photon. To 
a very large extend we can consider this state just a bound state of an electron, even 
if technically it should be considered a resonance because it couples to channels in the 
continuum at the states gets a width. 

The width of the resonance can be calculated from Eq. (1791) . To simplify the formulation 
let us assume that we have the bound channel a and the open channel 1 and neglect v n 
versus V\ a and v aa . Eq. ( 120]) gives now: 



1 



1 — v aa G a — V la GiG a 



v\ a G a Vi a 

Via vf a Gi + V a 

V n „ 



1 - v aa Ga - dG a /dE v aa {E - E R ) - vl a G x G a 

E=E K 



15) 

16) 



E=E R 



where, since for the energy E R we have approximately a bound state of a, according to [19] 
we will have v aa G a = 1 and hence 



■1 



9a 



dGJdE E £n vi„.Gf« a 



Vaa{-dGa/dE) 



E~E K + lTi/2 



^7) 



E=E R 



which leads to 



9a 



1 



-dGa/dE 



£i = ImivlG^a) 
2 



8Gg 

1 dE 



E=E R 

Considering that from Eq. ([9]) we find 

Im G\ = —An 2 niki . 
From Eqs. (J75JI . (BOP and (|ggj) we find 



E—E-b 



^ 1 a 1 i V la9 a G a A 2 ; 9l V aaG a 

= Air /iifci 

Vaa V a a 



4:Ti 2 fiik 



4-7T (iifagl 



19) 



(90) 



where in the last step of Eq. (190]) we have used Eq. ({84]) and v aa G a = 1 for the bound 
state. Eq. ( 190]) is the same as Eq. ( 166]) upon establishing the equivalence of the couplings 
in Field Theory and Quantum Mechanics of Eq. (160]) for i — 1, 

g\ = 32 7 r 3 //i^(2M B i)- 1 ^ . (91) 
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Another way to proceed, which stresses the picture that we have, is a straightforward 
derivation of the probability to decay into the channel j of the bound state a. The standard 
procedure [34] starts from the S matrix 



cc 



S = l-i Vidt + --- 

J — oo 

/oo 
e- iE ^e iEat (j\V\a)dt 
-oo 

= -i2irS(Ej - E a ) J d 3 p J d 3 p'(j\p)(p\V\p')(p'\a) 

-i2nS(E j - E a )v ja I d 3 p'(p'\a) = -i2irS{E j - E a ) Vj J a (92) 

Jv'<A 



where we have used J <A d 3 p(j\p) = 1 for the asymptotic free \j) state (\Pj)) and must 
remember to sum over \j) in \S\ 2 by means of d 3 pj to be consistent with our normalization. 
With the standard procedure to deal with (S(E — Ej)) 2 we find 

T -S d3p ' 2 ^- E - )V -i^r (93) 

where we have divided by the norm squared of the wave function of the bound state a. 
From [19] we know that for the bound state 

m^Pa)= d 3 p{ _ l _ ) 2 J2 V <d1>l Yl Vam ^rn (94) 

Jp<A ^a lvl a P I ^•H'a , 

which, upon the approximations used that |f aa | 3> |%| (i, j ^ a), reads as 

dG 

(lpa\lpa) = ^7T VaatpaVaatpl (95) 

ahi E=E a 

\ 1 _1 

and using the fact that g 2 = (—dG a /dE) E=Ea we find 

2 
■n o 2 9 a V ja V ja o 2 2 fnr\ 

T = Sir HjPj-*^ — — = 57T fijPjgj (96) 

where in the last step we used Eq. ( 1841) . Thus, the result is the same as in Eq. ( 190~1) . 

The derivation has served to see that this result holds in the case that a is a domi- 
nant channel with moderate decay into open ones, providing this intuitive picture for the 
resonant states in coupled channels. 
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3 Application to the two A(1405) states 

3.1 The A(1405) and A(1670) resonances in the Chiral Unitary 
Model 

We apply the formalism explained in the former sections to study the wave functions 
of the two A(1405) states generated in the chiral unitary approach [14,15,23]. For the 
description of the A(1405) states, we use the chiral unitary approach [23,24]. In this model 
there are two A (1405) states dynamically generated in the coupled channels of meson- 
baryon scattering, KN, 7r£, T]h and KE. The most interesting thing in this model is 
that two poles exist around the A (1405) energy region at z\ = (1390, — i66) MeV and 
z 2 = (1426, — il6) MeV [23]. These two poles are also found in all the works on the chiral 
unitary approach [35-37] that followed Ref. [23], even when higher order terms in the 
chiral Lagrangians are considered. Experimental support for these states has been shown 
in [38,39]. The electromagnetic mean squared radii of A(1405) are calculated in Ref. [40] 
and are shown to be much larger than that of ground state baryons. The A(1670) was first 
reported as a dynamically generated resonance in Refs. [15,24] and has been corroborated 
in following works [14,17]. 

3.2 Wave function in coordinate space 

We study here the wave functions in / = of the KN, 7rS, r]A, KE channels. Since the 
A (1405) is observed in the 7rS spectra in experiment, we consider that the 7rS channel is 
the scattering state that we have called channel 1. 

We show in Fig. [3] the wave function in coordinate space at the pole energies of the 
two A(1405), together with that of the A(1670). From Fig. El we find that the KN 
components dominate at 1426 MeV while the 7rS components are dominant at 1390 MeV. 
This is consistent with the findings of [14,15,23,35-37] that the pole at higher energies 
couples most strongly to KN while the one at lower energies couples mostly to 7rS. 

The curves correspond to these different energies, which are the energies where we 
find the poles of the A(1390), A(1426) and A(1670). We can see how the wave function 
concentrates close to the origin, and both for bound channels, as well as for open channels, 
fades away rapidly beyond 2 fm, providing a spatial picture of the distribution of the 
particles of the different channels building up the resonances. 

In the case of the A(1670) we can see in the last line of Fig. [3] that the dominant 
component is the KE bound state. This resonance would very approximately qualify as a 
KE bound state, as also suggested in Ref. [24] based on the large couplings of the resonance 
to that state. 
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Figure 3: Wave functions in coordinate space. Solid lines and dashed lines show the results 
obtained at E = 1426 MeV and at E = 1390 MeV for the A(1405) resonances, and dotted 
lines show that at E — 1680 MeV for the A(1670) state. The figures of the first line contain 
the full wave function of 7rS (initial plus scattered), while the second line shows only the 
scattered wave of the 7rS state 
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4 Response function and form factors 

As a further application of the formalism, let us compute the response function of one 
resonance state to an external scalar source. Let us take one of the components (one 
channel), the final expression will contain the sum of the partial response functions weighted 
by its coupling to the resonance squared, as we shall see. Let us assume the channel to be 
a bound state for the moment, and let us couple the external scalar source, with strength 
unity, to one of the particles, particle 1 in Fig. HI The S matrix for this diagram is given 



\ 
k' ^ 



\ 



/ 



/i 



/ 



Xi 



x 2 



A 



by 



Figure 4: Coupling on external source to a molecular component 



S 



d A Xi 



\A^ 



- ip ° x °M Xl ) - 



VI 



JEi. 



=e^0!« 



-ikx\ 



'A 



VzJkV V^^v 



e ik ' xi (-iti 



(97) 



where V is the volume of a box where we normalize to unity our wave functions. In Eq. (|9 
t\ = 1 for our scalar source. We can multiply by 



d 3 X 2 4>2{X2)4>2(X2) = 1 



(98) 



and perform the x\ integration which provides the 6 of conservation of energy. Furthermore 

we can now write 

(99) 



h{x 1 )<p 2 {x 2 ) 



V 



e iKm - R 6(x) 



where R, x are the CM and relative coordinates given by 

vf\\X\ + m 2 X2 



R 



x 



m\ + m 2 
x 2 - Xi 



(100) 
(101) 
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and K m , <p are the total momentum of the molecule and its relative wave function. 
After performing the d 3 x\d 3 x 2 — > d 3 Rd 3 x integrations we obtain 

S = -iti±i 1— J— 1- J_ (2vr)V(fe + K m -k'- K'JF 
where F(q) is the form factor given by 



iqx 



m 2 
m\ + m 2 


(k- 


(102) 
(103) 


o _ Jd 
n — Pi 


- K°- 


The 



assuming the wave function real for the bound state, and K^ — K n 
scalar source can couple to particle 2 and we would have to sum this contribution, which 
is trivially obtained by exchanging the 1 and 2 indices in the former expressions. 
It is useful to go to momentum space to evaluate Eq. (11031) . We have 

<K*) = W> = /(^ e ^<2#> (104) 

(pm = v 6 ^~ p) [ d 3 k(k\^} . (los) 

E-Uxip) -uj 2 {p) J k<A 

Hence, although we can take 4>(x) real for a bound components (we will see latter on the 
generalization to open channels), formally we can write 

F(,) = /^(^»e- 

d'p 6(A-p) 



d x I - — rTT^e 



(2tt)3/2 E - u^p) - u 2 (p) 

d 3 p' c - ip >. x e(A-p') c _ iq . x 



(2tt)3/ 2 E - u^p') - u 2 {p>) 

d * p B{A-p)e{A-\p-g\) 

(E - wi(p) - cu 2 (p))(E - wi(p - q) - u 2 (p - q)) 

up to a normalization, easily restored demanding that F(q = 0) = 1. 

Let us now compare this result with what one would obtain in a field theoretical ap- 
proach in which the scalar source couples to the components of the resonance. Diagram- 
matically the process is depicted in Fig. The response function from interaction with 
particle 1 will be (taking P = 0, the total three momentum of the resonance) 



d'p 9{A-p) 0(A-|p-^g|) 1 

(27r) 4 p° — U)i{p) + iep° — q° — iO\{p — q) + ie P° — p° — u 2 {p) + ie 

which, upon, contour integration of the p° variable on the upper half circle gives 

d 3 P fl(A-p) g(A-lp-^^gl) 

(2vr)3 po _ W2 ( p ) _ Wl ( p ) + ieP o _ q o _ Wl ( p -q)- u 2 (p) + it 
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(107) 
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Figure 5: 



with q° the excitation energy carried by the external source in this reference frame, such 
that P°-q° = E R (q), q° = M R - E R (q) = -q 2 /2( mi + m 2 ). 

Even if Eq. (11061) and Eq. (11081) do not look the same, one can see that in the nonrela- 
tivistic limit one has 

p x = miii; p 2 = m 2 x 2 ; P = MR = p x +p 2 

H= ; , P = ^ x = -77P2 ~ 77P1 , (109) 

Wj\ + m 2 M M 

and then the field theoretical approach with the loop function simply provides the form 
factor that one obtains in the wave function approach. The argument of the second 9 
function in Eqs. (11071) and (11081) corresponds to the relative momentum for particles 1 and 
2, having momenta p — q and —p, as shown as Fig. El which is — p H ^r—q- 

ifl 2 1* lit 2 

One subtlety is worth mentioning at this point. In the case of open channels in the wave 
function method, the evaluation of the form factor with the used f <f)(x)<f)*(x) exp(— iq ■ 
x)d 3 x would lead to the expression of Eq. f |108[) with a — ie instead of ie in the second 
factor. The field theoretical approach, which provides an appropriate formalism for these 
processes, keeps the +ie in the two propagators. This means that in the case of interaction 
of an external source with the open channels of a resonance the response function does 
not involve the ordinary form factor involving <jxfr* but something else, which in the case 
of bound channels is the ordinary form factor. The situation can also be interpreted from 
the quantum mechanical side, since in the case of decay into an open channel, the 0* 
conjugate of the wave function of the final state is an outgoing solution of the Schrodinger 
equation with the potential involved (complex to account for inelastic channels) and not of 
the complex conjugate of the Schrodinger equation [41]. In the case that we have several 
channels building up a resonance, the response function is given by 

R(q)= Yl 9iGi(E,q) gi (110) 

channels 
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with g.i the coupling of the resonance to each channel and 



G, 



d 3 p 



6(A-p)6(A- \p 



llt'2, 



m li +m 2i 



Q)\) 



(2tt) 3 / 2 (e - u;i(p) - w 2 (p) + ie){^E* + q 2 - u t (p - q) - co 2 (p) + it) 



111^ 



In the response function we would have to sum over the second particle too. 
We show results in Figs. M to M 
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Figure 6: Contribution of the meson part of each channel to the response function for the 

two A(1405) states. 



In Fig. Owe show the contribution to the response function from the meson component 
of each channel for two different energies, which correspond to the two A (1405) resonances 
of our approach. We see a gradual fall down of the response function with q, approaching 
zero at momentum of the order of 1 GeV. In Fig. [7| we show the corresponds contributions 
from the baryon part of the resonances. Finally in Figs. [SJ [9] we show the sum of the 
contributions of the different channels from the meson and baryon parts respectively. One 
observes differences in the behavior of the scalar response function for the two resonances, 
with a faster fall down for the A (1420) case. Although the scalar form factor is different 
than the charge form factor, it is interesting to note that there are similarities of our results 
with the charge form factors evaluated in [40], where the one of the A(1420) has also a 
faster fall down with q than the one of the A(1390). 
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Figure 7: Contribution of the baryon part of each channel to the response function for the 
two A(1405) states. 
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Figure 8: Sum of all contribution from the meson part of the two resonances. 
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Figure 9: Sum of all contribution from the baryon part of the two resonances. 



5 Summary 

In this paper we have developed a formalism to deal with coupled channels in a unitary 
approach by paying a special attention to constructing the wave functions in the different 
channels in the case that there are resonances dynamically generated. The paper general- 
izes what was found before for only bound coupled channels. Here we have bound and open 
channels and the formalism is subtly different, since contrary to the case of bound states, 
where only discrete energies are allowed, here we have a continuous energy variable. Many 
of the results obtained for bound states do not hold for the resonance states. One of the 
things we do is to identify the meaning of a resonance in the coupled channel approach, and 
it emerges as an approximate bound state of a coupled channel which can decay into the 
open ones. The formalism developed is easy, practical and useful. A separable potential in 
coordinate space is chosen which leads to an on shell factorization of the Bethe Salpeter 
equations (Lippmann Schwinger in the nonrelativistic form), which allows to convert the 
coupled channel integral equations into trivial algebraic equations. The wave functions in 
momentum space are then found as trivial analytic functions, from where the wave func- 
tions in coordinate space can be easily evaluated. The couplings of the resonance to the 
different channels are related to the wave function at the origin and interesting relation- 
ships between these couplings are obtained. We also study the issue of couplings of the 
resonances to states outside the space of the building channels and justify results used 
before in the Literature, setting the limits for their application. Similarly, we also face 
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the issue of final state interaction within the coupled channel formalism and find again a 
justification for results used in the Literature, setting again the limits of applicability. 

As an application of the formalism, we tackle the problem of the two A(1405) and the 
A(1670) states dynamically generated in the chiral unitary approach from the 7rS, KN, rjA, 
and ICE interaction. We evaluate the wave functions in coordinate space for the first time, 
giving an intuitive idea of the wave functions and the spatial distribution of the particles 
of the different channels. We have also evaluated the response function of the resonances 
to an external scalar source, together with the form factors involved. 

We envisage practical applications of the formalism in any problem where one has to 
deal with the interaction of external sources with the dynamically generated resonances. 
In particular in the building up of multihadron states using iteratively the Fixed Center 
Approximation to the Faddeev equations, where other methods would turn out technically 
prohibitive. The resulting wave functions can also be used to evaluate static properties of 
the resonances and different form factors, etc, which might be studied experimentally in 
the future. 
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